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Abstract 



Using the effective potential approach for composite operators, we have for- 
mulated a general method how to calculate the truly nonperturbative vacuum 
energy density in the axial gauge QCD quantum models of its ground state. 
It is defined as integrated out the truly nonperturbative part of the full gluon 
propagator over the deep infrared region (soft momentum region). The non- 
trivial minimization procedure makes it possible to determine the value of 
the soft cutoff in terms of the corresponding nonperturbative scale parameter 
which is inevitably presented in any nonperturbative model for the full gluon 
propagator. If the chosen Ansatz for the full gluon propagator is a realistic 
one, then our method uniquely determines the truly vacuum energy density, 
which is always finite, automatically negative and it has no imaginary part 
(stable vacuum). We illustrate it by considering the Abelian Higgs model of 
dual QCD ground state. We have explicitly shown that the vacuum of this 
model without string contributions is unstable against quantum corrections. 

PACS numbers: ll.15.Tk, 12.38.Lg 
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I. INTRODUCTION 



The nonperturbative QCD vacuum is a very complicated medium and its dynamical 
and topological complexity [1-3] means that its structure can be organized at various levels 
(classical, quantum) and it can contain many different components and ingredients which 
contribute to the vacuum energy density (VED), one of the main characterics of the QCD 
ground state. Many models of the QCD vacuum involve some extra classical color field 
configurations such as randomly oriented domains of constant color magnetic fields, back- 
ground gauge fields, averaged over spin and color, stochastic colored background fields, etc. 
(see Refs. [1,4] and references therein). The most elaborated classical models are random 
and interacting instanton liquid models (RILM and IILM, respectively) of the QCD vacuum 
[5]. These models are based on the existence of the topologically nontrivial instanton- type 
fluctuations of gluon fields, which are nonperturbative solutions to the classical equations of 
motion in Euclidean space (see Ref. [5] and references therein). 

Here we are going to discuss the quantum part of VED which is determined by the 
effective potential approach for composite operators [6,7] (see also Ref. [8]). It allows us 
to investigate the nonperturbative QCD vacuum, in particular Yang-Mills (YM) one, by 
substituting some physically well-justified Ansatz for the full gluon propagator since the 
exact solutions are not known. In the absence of external sources the effective potential is 
nothing but VED which is given in the form of the loop expansion where the number of the 
vacuum loops (consisting in general of the confining quarks and nonperturbative gluons) is 
equal to the power of the Plank constant, h. 

Let us remind the reader that the full dynamical information of any quantum gauge field 
theory such as QCD is contained in the corresponding quantum equations of motion, the 
so-called Schwinger-Dyson (SD) equations for lower (propagators) and higher (vertices and 
kernels) Green's functions [9,10]. These equations should be also complemented by the cor- 
responding Slavnov- Taylor (ST) identites [9-12] which in general relate the above mentioned 
lower and higher Green's functions to each other. These identities are consequences of the 
exact gauge invariance and therefore "are exact constaints on any solution to QCD" [9]. 
Precisely this system of equations can serve as an adequate and effective tool for the nonper- 
turbative approach to QCD. Among the above-mentioned Green's functions, the two-point 
Green's function describing the full gluon propagator (see section II below) has a central 
place [9-13]. In particular, the solutions to the above-mentioned SD equation for the full 
gluon propagator, are supposed to reflect the quantum structure of the QCD ground state. 
It is a highly nonlinear integral equation containing many different propagators, vertices 
and kernels [9-13]. For this reason it may have many different exact solutions with different 
asymptotics in the deep infrared (IR) limit (the ultraviolet (UV) asymptotics because of 
asymptotic freedom are apparently uniquely determined), describing thus many different 
types of quantum excitations of gluon field configurations in the QCD vacuum. Evidently, 
there is no hope for an exact solutions as well as not all of them can reflect the real structure 
of the QCD vacuum. Let us emphasize now that any deviation in the behavior of the full 
gluon propagator in the IR domain from the free one automatically assumes its dependence 
on a scale parameter (at least one) resposible for nonperturbative dynamcs in the quntum 
model under consideration, say, An p. This is very similar to asymtotic freedom which 
requires asymptotic scale parameter associated with the nontrivial perturbative dynamics 
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(scale violation). However, to calculate the truly nonperturbative VED we need not the IR 
part in the decomposition of the full gluon propagator, but rather its truly nonperturbative 
part which wanishes when the above-mentioned nonperturbative scale parameter goes to 
zero, i.e., when the perturbative phase survives only in the corresponding decomposition of 
the full gluon propagator (see next section below). 

It is well known, however, that VED is badly divergent in quantum field theory, in 
particular QCD (see, for example, the discussion given by Shifman in Ref. [1]). The main 
problem thus is how to extract the truly nonperturbative VED which is relevant for the 
QCD vacuum quantum model under consideration. It should be finite, negative and it 
should have no imaginary part (stable vacuum). Why is it so important to calculate it from 
first principles? As was emphasized above, this quantity is important in its own right being 
nothing but the bag constant (the so-called bag pressure) apart from the sign, by definition 
[14]. Through the trace anomaly relation [15] it assists in the correct estimating such an 
important phenomenological nonperturbative parameter as the gluon condensate introduced 
in the QCD sum rules approach to resonance physics [16]. Furthermore, it assists in the 
resolution of the U(l) problem [17] via the Witten-Veneziano (WV) formula for the mass of 
rf meson [18]. The problem is that the topological susceptibility needed for this purpose [16- 
19] is determined by the two point correlation function from which perturbative contribution 
is already subtracted by definition [18-22]. The same is valid for the above-mentioned bag 
constant which is much more general quantity than the string tension since it is relevant for 
light quarks as well. Thus to correctly calculate the truly nonperturbative VED means to 
correctly understand the structure of the QCD vacuum in different models. 

We have already formulate a method how to calculate the truly nonperturbative YM 
VED in the covariant gauge QCD [23]. The main purpose of this paper (section II) is to 
formulate precisely a general method how to correctly calculate the truly nonperturbative 
quantum part of YM VED in the axial gauge QCD. In sections III and IV we illustrate 
it by considering the Abelian Higgs model [24] of the dual QCD [25] ground state. We 
will explicitly show that the vacuum of this model without string contributions is unstable 
against quantum corrections. In section V we summarize our results. 

II. THE TRULY NONPERTURBATIVE VACUUM ENERGY DENSITY 

In this section we are going to analytically formulate a general method of calculation 
of the quantum part of the truly nonperturbative YM VED in the axial gauge QCD. Let 
us start from the nonperturbative gluon part of VED which to-leading order (log-loop level 
~ h)^\ is given by the effective potential for composite operators [6] as follows 

V{D) = \S w^ Tr[HD ° lD) ~ {D ° lD) + 1} ' (2 - 1} 

where D(q) is the full gluon propagator (see below) and D (q) is its free (perturbative) coun- 
terpart. Here and below the traces over space-time and color group indices are understood. 



1 Next-to- leading and higher terms (two and more vacuum loops) are suppressed by one order of 
magnitude in powers of fi at least and are left for consideration elsewhere. 
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The effective potential is normalized as V(D ) = 0, i.e., the free perturbative vacuum is 
normalized to zero. 

A general parametrization of the gauge boson propagator in the axial gauge of dual QCD 
is [24-26] (here and below we use notations and definitions of Refs. [24,26]) 



1 



(q-n) 



;T, u (n)G(-q 2 ) + V(g,n)F(-g 2 ), 



(2.2) 



where 



Luv(q,n) = 5„„ - 



q^n v + q u n^ q^q u 



+ 



(q -n) (q- n) 



(2.3) 



with an arbitrary constant unit vector n M , n 2 ^ = 1. The exact coefficient functions G(—q 2 ) 
and F(—q 2 ) characterize the vacuum of the theory under consideraion. Their free perturba- 
tive counterparts are 



?PT, 



~q 2 ) 



G PT (~q 2 



o. 



(2.4) 



Thus the free perturbative gluon propagator is 

1 



D%(q,n) 



-q 2 ) 



L^u(q,n) 



while its inverse is 



"" ~ q 2 , 



(2.5) 



(2.6) 



Using futher Eqs. (2.2) and (2.6), one obtains 



[DL\ mA^n) = (-q 2 )F(-q 2 ) + G(-q 2 ). 



(2.7) 



In order to evaluate the effective potential (2.1) we use the well-known expression, 



Tr H D o lD ) = 8 x IndetiD^D) = 8 x 4 In {-q 2 )F{-q 2 ) + G{~q 2 



(2.8) 



It becomes zero (in accordance with the above mentioned normalization condition) when 
the full gluon form factors are replaced by their free counterparts (see Eqs. (2.4)). Going 
over to four {n = 4) dimensional Euclidean space in Eq. (2.1), on account of (2.8), and 
evaluating some numerical factors, one obtains (e 9 = V(D)) 

e 9 = ~^J dq 2 q 2 [in (q 2 F(q 2 ) + G(q 2 )) - (q 2 F(q 2 ) + G(q 2 )) + l] . (2.9) 

Let us now introduce the following decomposition of the exact coefficient functions G(q 2 ) 
and F(q 2 ) (Euclidean metrics) 
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Ftf) = F»"( t t t ) + F"tf), 

G(q 2 ) = G NP (q 2 ) + G PT (q 2 ), (2.10) 
where the truly nonperturbative quantities F NP (q 2 ) and G NP (q 2 ) are defined as follows: 



F NP ( q 2 ,A NP ) 

G NP (q 2 ,A NP ) 



F(q 2 ,A NP )-F(q 2 ,A NP = 0), 
G(q 2 ,A NP ) ~G(q 2 , A NP = 0), 



(2.11) 



which explains the difference between the truly nonperturbative parts and the full gluon 
form factors which are nonperturbative themselves. Let us note, that the perturbative 
parts F PT (q 2 ) and G PT (q 2 ) may, in general, contain renormgroup log improvmenets due 
to asymptotic freedom. Without these improvments their free perturbative counterparts 
are given in Eqs. (2.4). Substituting these relations into Eq. (2.9) and doing some trivial 
rearangment, one obtains 



1 



e n = — 



7T Z 



dq 2 q 2 [in (l + q 2 F NP (q 2 ) + G NP (q 2 )) - (q 2 F NP (q 2 ) + G NP (q' 



+ IpT, 

(2.12) 



where we introduce the following notation 



i PT 



= o dq 2 q 



7T" 



In 1 



1 - q F (q 



G 



PT 



1 + q 2 F NP {q 2 ) + G NP {q 2 ) 




1 - q 2 F PT (q 2 ) - G 



PT 



(S 2 )) 

(2.13) 



as containing contribution which is mainly determined by the perturbative part. However, 
this is not the whole story yet. We must now to introduse the soft cutoff in order to 
separate the deep IR region where the truly nonperturbative contributions become dominant 
(obviously they can not be valid in the whole energy- momentum range). So the expression 
(2.12) becomes 



~ [ l dq 2 q 2 [in (l + q 2 F NP {q 2 ) + G NP {q 2 )) - (q 2 F NP {q 2 ) + G NP {q 2 )) 



+ IpT + IpTi 

(2.14) 



where the explicit formula for I PP (which is obvious) is not important. The contribution 
over the perturbative region I PT as well as I PP should be subtracted by introducing the 
corresponding counter terms into the effective potential, which is equivalent to define the 
truly nonperturbative VED as e" p = e g — I PT — I PT - Thus one finally obtains 



[ l dq 2 q 2 [in (l + q 2 F NP {q 2 ) + G NP {q 2 )) - (q 2 F NP (q 2 ) + G NP (q 2 ))} . (2.15) 



This a general formula which can be applied to any model of the axial gauge QCD ground 
state based on the corresponding Ansatz for the full gluon propagator. So Eq. (2.15) is our 
definition of the truly nonperturbative VED as integrated out the truly nonperturbative part 
of the full gluon propagator over the deep IR region, soft momentum region, < q 2 < q 2 . 
How to determine q^l By the corresponding minimization procedure, of course (see below). 
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A. 

From this point it is convenient to factorize the dependence on a scale in the nonpertur- 
bative VED (2.15). As was already emphasized above, the full gluon form factors always 
contain at least one scale parameter responsible for the nonperturbative dynamics in the 
model under consideration, A^rp- Within our general method we are considering it as free 
one, i.e., as "running" (when it formally goes to zero, only perturbative phase survives in 
the model under consideration) and its numerical value (if any) will be used only at final 
stage in order to numerically evaluate the corresponding truly nonperturbative VED (if any). 
We can introduce dimensionless variables and parameters by using completely extra scale 
(which is aways fixed in comparison with Atvp), for example flavorless QCD asymptotic scale 
parameter Aym as follows: 

a 2 a 2 A 2 

Z = T^> ^o = -r^, b =T^- ( 2 - 16 ) 

IV YM IV YM IV YM 

Here zq is the corresponding dimensionless soft cutoff while the parameter b has a very clear 
physical meaning. It measures the ratio between nonperturbative dynamics, symbolized by 
A^p and nontrivial perturbative dynamics (violation of scale, asymptotic freedom) symbol- 
ized by Ay M . When it is zero only perturbative phase remains in the quantum model under 
consideration. In this case, the gluon form factors obviously become a functions of z and 
6, i.e., F NP (q 2 ) = F NP (z,b) and G NP (q 2 ) = G NP (z,b), so the truly nonperturbative VED 
(2.15) is (e n / = e^(z ,b)) 

n g (z ,b) = -^e^(z ,b), (2.17) 

ly YM 

where for futher aims we introduce the gluon effective potential at a fixed scale Ay M , [23,27] 

fl g = Q g (z , b) = ^ [J dz z [(zF NP (z, b) + G NP (z, b)) - In (l + zF NP (z, b) + G NP (z, b 

(2.18) 

Precisely this expression allows us to investigate the dynamical structure of the YM vacuum 
free of scale dependence complications as it has been already factorized in Eq. (2.17). It 
depends only on zo and b and the minimization procedure can be done now with respect to b, 
dQ g (z , b)/db — (usually after integrated out in Eq. (2.18)) in order to find self-consistent 
relation between zo and b, which means to find go as a function of A^p. Let us note in 
advance that all final numerical results will always depend only on A^p as it should be for 
the nonperturbative part of VED. Obviously, the minimization with respect to z leads to 
trivial zero. In principle, through the relation A YM = q^ZQ 2 , it is possible to fix the soft 
cutoff q itself, but this is not the case indeed since then z can not be varied. 

B. 

On the other hand, the scale dependence can be factorized as follows: 
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A 2 ' 
NP 



Z 



A 2 ' 

NP 



(2.19) 



i.e., 6=1. For simplicity (but not loosing generality) we use the same notations for the 
dimensionless set of variables and parameters as in Eq. (2.16). In this case, the gluon form 
factors obviously becomes the function of z only and the truly nonperturbative VED (2.15) 
becomes 

€ 7( z °) = ^9oV £ dz z [(zF NP (z) + G NP (z)) - In (l + zF NP (z) + G NP (z))} . (2.20) 

Evidently, to fix the scale now is possible in the two different ways. In principle, we can fix 
A NP itself, i.e., introducing 

Q g (z ) = xrt^o) = ^ jf dz z [(zF NP (z) + G NP (z)) - In (l + zF NP (z) + G NP (z)) 

(2.21) 

However, the minimization procedure again leads to the trivial zero, which shows that this 
scale can not be fixed. 

In contrast to the previous case, let us fix the soft cutoff itself, i.e., setting [23,28] 



= Jf7^) = ^ [° dz z [(zF NP (z) + G NP (z)) - In (l + zF NP (z) + G NP (z))'_ 

(2.22 



The minimization procedure with respect to z is nontrivial now. Indeed, dQ g (z )/dz = 0, 
yields the following "stationary" condition 



j*° dz z [(zF NP {z) + G NP {z)) - In (l + zF NP (z) + G NP (z)) 
(z F NP (z ) + G NP (z )) - In (l + z F NP (z ) + G NP (z 



(2.23) 

which solutions (if any) allows one to find qo as a function of A^p. On account of this 
"stationary" condition, the effective potential (2.22) itself becomes simpler for numerical 
calculations, namely 

n,(tf) = ^ [{*#F» p {z*) + G NP (tf)) - In (l + 4F NP {4) + G NP (z«))] , (2.24) 

where Zq is a solution (if any) of the "stationary" condition (2.23) and corresponds to the 
minimum(s) (if any) of the effective potential (2.22). In the next sections we will illustrate 
how this method works. 



III. ABELIAN HIGGS MODEL 

Let us now consider some special model of the dual QCD [25] ground state. In the 
dual Abelian Higgs theory which confines electric charges the coefficient functions F(q 2 ) 
and G(q 2 ) are [24] (Euclidean metrics) 
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M B D*(q 



q 2 + M B \ q 2 + M% J 
G(a 2 )~ Ml (l M 2 <fD ^ q2) \ (3 1) 

where M B is the mass of the dual gauge boson and D T '(q 2 ) represents the string contribu- 
tion into the gauge boson propagator. The mass scale paprameter M B is the scale responsible 
for nonperturbative dynamics in this model (in our notations Ajyp = M B ). When it formally 
goes to zero, then one recovers the free perturbative expressions indeed, (2.4). Removing 
the string contributions from these relations we get 

1 M 2 

jpno—str. I „2\ /~mo—str.(„2\ lv± B fo Q\ 

F (g) = ^TM|' G (g) = - ?TM|' (3 ' 2) 

i.e., even in this case these quantities remain nonperturbative. The truly nonperturbative 
expressions (2.11) now become 

F NP (n 2, = M 2 B { _ M 2 B q 2 D\q 2 ) \ 

^ ] ~ Q 2 (q 2 + M 2 ){ q 2 + M 2 J' 

G {q) -~^MlV~ q 2 + M 2 B J' (3 ' 3) 
while with no-string contributions they are 

M 2 M 2 

Tpno—str.f J2\ 1V± B r^no—str. ( „2\ 1V± B (q a \ 

FnP {q) ~'q 2 (q 2 + M 2 B y ° NP (q ) -~W+Ml' ( ^ 

Both expressions (3.3) and (3.4) are truly nonperturbative indeed, since they become zero 
in the perturbative limit (M B — > 0), when only perturbative phase remains. From these 
relations also follows 

r NP ( 2, 2 F NP ( 2x_ Ml ( M 2 B q 2 D^q 2 ) \ 
G (q)-qF (, ) - --^—^l - ?2 + m| J, 

M 2 

G n N V tr \l 2 ) = q 2 F n N °f str \q 2 ) = --TTT72> ( 3 - 5 ) 

q -\- ivi B 

so the truly nonperturbative vacuum energy density (2.15) will depend only on one function, 
say, G NP (q 2 ) (see next section). 

Although the expression (2.2), on account of (3.1), for the gluon propagator is exact, 
nevertheless it contains an unknown function D s (g 2 ) which is the intermadiate string state 
contribution into the gauge boson propagator [24]. It can be considered as a glueball state 
with the photon quantum numbers 1 _ . The bahavior of this function D s (q 2 ) in the IR 
region (q 2 — > 0) can be estimated as follows [24]: 

^ 2 H^ + ..., (3.6) 
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where C is a dimensionless parameter and M gl is the mass of the lowest 1~ glueball state. 
The dots denote the contributions of heavier states. Thus, according to Eqs. (3.1) and (3.6), 
the coefficient functions in the IR limit behave like 

F ^) = W B + ^ + 0{q% G(q 2 ) = -l + 0(q 2 ), g 2 -0. (3.7) 

At the same time according to Eqs. (3.3), (3.5) and (3.6) their truly nonperturbative coun- 
terparts behave like G NP (q 2 ) = q 2 F NP (q 2 ) = — 1 + 0(q 2 ), q 2 — > 0, i.e., in the same way as 
G{q 2 ) in Eq. (3.7). 



IV. VACUUM STRUCTURE IN THE ABELIAN HIGGS MODEL 

Let us calculate the truly nonperturbative VED in the Abelian Higgs model described 
in the preceding section. It is instructive to start from the case when there are no string 
contributions into the structure functions F(q 2 ) and G(q 2 ). Then their truly nonperturbative 
parts are given in Eqs. (3.4). It is convenient to factorize the scale dependence of VED by 
introducing dimensionless variables and parameters in accordance with B-scheme (2.19) with 
A NP = M B - In this case, the gluon form factors (structure functions) obviously becomes the 
functions of z only, q 2 F NP {q 2 ) = G NP {q 2 ) = -(1/1 + z). The truly nonperturbative VED 
(2.22) becomes 



- I i r z o 

^o) = -iC( 2 o) = -- 2 *» 2 / dz z 

yg 71 J 

Easily integrating Eq. (4.1), one obtains 
1 



1 + z 



+ In 



-1+z 
l + z 



(4.1) 



n g (z ) 



'2tt 2V 



2^ -41n(l + ^ )+ln(l±^)+^ 2 h/ ' '" ' 



1 - z 



1+Zq 



(4.2) 



From this expression it is almost obvious that the effective potential will have imaginary 
part at any finite value of the soft cutoff, which is a direct manifestation of the vacuum 
instability [29]. Asymptotics of the effective potential (4.2) to-leading order are 



n g (z ) 



2 ^0 



2tt 2 



ln(-l), 



n g (z ) 



— z 2 lnz . 



7T Z 



Let us remind, that z — > oo is the perturbative limit (Mb 
fixed. The "stationary" condition (2.23) now is 



(4.3) 

0) when the soft cutoff q is 



41n(l + z ) - In 



l + z 



((1 + z ) 2 + 2) . 



. 1 - z J 1 + z 

It has only trivial solution z = 0, so the "stationary" state does not exist in this model 



(4.4) 
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V. CONCLUSIONS 



In summary, we have formulated a general method how to calculate the truly nonpertur- 
bative VED in the axial gauge QCD quantum models of its ground state using the effective 
potential approach for composite operators. It is defined as integrated out the truly nonper- 
turbative part of the full gluon propagator over the deep IR region (soft momentum region). 
The nontrivial minimization procedure which can be done only by the two different ways 
(leading however to the same numerical value (if any) of VED) makes it possible to determine 
the value of the soft cutoff in terms of the corresponding nonperturbative scale parameter 
which is inevitably presented in any nonperturbative model for the full gluon propagator. If 
the chosen Ansatz for the full gluon propagator is a realistic one, then our method uniquely 
determines the truly nonperturbative VED, which is always finite, automatically negative 
and it has no imaginary part (see, for example our previous publications [23,28]). Here 
we illustrate it by considering the Abelian Higgs model of the dual QCD ground state. 
The quantum part of VED (4.2) always contains imaginary part. Thus, the vacuum of the 
Abelian Higgs model without string contributions is unstable indeed. Whether the string 
contributions can cure this fundamental problem or not is beyond the scope of this letter 
and is left for consideration elsewhere. The vacuum instability of the Abelian Higgs model 
without string contributions will be recovered, of course, within A-scheme (2.16) as well. 
Nothing depends on how one introduces the scale dependence by chosing different scale 
parameters. 

Comparing Eqs. (2.9) and (2.15), a prescription to obtain the relevant expression for 
the truly nonperturbative VED can be derived. Indeed, for this purpose in Eq. (2.9) the 
replacement q 2 F(q 2 )+G(q 2 ) -> l+q 2 F NP (q 2 ) +G NP (q 2 ) should be done. Also the soft cutoff 
q 2 on the upper limit should be introduced. Now it looks like the UV cutoff, but nevertheless 
let us underline once more that it separates the deep IR region from the perturbative one, 
which includes the IM region as well. It can not be arbitrary large as the UV cutoff is, 
by definition. As far as one chooses Ansatz for the full gluon propagator, the separation 
"NP versus PT" in Eq. (2.10) is exact because of the definition (2.11). The separation "soft 
versus hard" momenta is also exact because of the above-mentioned minimization procedure. 
Thus the proposed determination of the truly nonperturbative VED is uniquely defined. It 
is possible to minimize the effective potential at a fixed scale (2.17) with respect to the 
physically meaningful parameter. When it is zero, the perturbative phase only survives in 
all quantum models of the QCD ground state. Equivalently, we can minimize the auxiliarly 
effective potential (2.22) as a function of the soft cutoff itself. As was underlined above, 
both methods lead to the same numerical value for the truly nonperturbative VED. 

There is no general method to formulate in order to calculate the confining quark quan- 
tum contribution into the total VED since this contribution depends heavily on the particular 
solutions to the quark SD equation. If it is correctly calculated then it is of opposite sign 
to the nonperturbative gluon part and it is one order of magnitude less as well (see, for 
example our recent papers [19,23,28]). Concluding, let us note that the generalization of our 
method on different noncovariant gauges [26,30] is straightforward. Let us underline that 
our method is not a solution to the above-mentioned fundamental badly divergent problem 
of VED. However, it is a general one and can be applied to any nontrivial QCD quantum 
vacuum models in order to extract the finite part of the truly nonperturbative VED in a 
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self-consistent way. In particular, it can serve as a test of different axial gauge QCD quan- 
tum as well as classical vacuum models since our method provides an exact criterion for the 
separation "stable versus unstable vacua". 
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